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CN . 1 Introduction, General Formalism 

O . 

D , The research on the superfluidity of neutron matter can be traced back to 

£^ ' Migdal's observation that neutron stars are good candidates for being macro- 

scopic superfiuid systems Q. And, in fact, during more than two decades of 
neutron-star physics the presence of neutron and proton superfiuid phases has 
been invoked to explain the dynamical and thermal evolution of a neutron star. 
The most striking evidence is given by post-glitch timing observations J^|| , but 
also the cooling history is strongly influenced by the possible presence of super- 
fluid phases |5|5|] . On the theoretical side, the onset of superfluidity in neutron 
, matter or in the more general context of nuclear matter was investigated soon 

^NJ ' after the formulation of the Bardeen, Cooper, and Schrieffer (BCS) theory of 

superconductivity |^] and the pairing theory in atomic nuclei [^|^] . 

The peculiar feature of a nucleon system is that it is a strongly interacting 
Fermi system with a force which has a short-range repulsive component and a 
' long-range attractive one. The first question raised by scientists was whether or 

not the strongly repulsive core might prevent the formation of a superfiuid state. 
But it was indeed shown |)| that the BCS approach, based on the mechanism of 
Cooper pairs, can be successfully extended to nuclear matter and that superfiuid 
states could in fact exist for a wide class of nucleon-nucleon potentials |l(| . The 
second question is related to the fact that the superfiuid state of nuclear matter 
is a selfsustaining state in the sense that nucleons participating to the pairing 
coupling also screen the pairing itself. From this point of view one expects the 
5J] ■ strong correlations to play an important role in delimiting the magnitude of the 

pairing gap. Therefore it appears necessary to go beyond the pure BCS approach 
and properly add the effects of the medium polarization as self-energy and vertex 
corrections. 

Since the neutron star "laboratory" can only provide an indirect evidence of 
the nucleon pairing in infinite matter and its relation with the pairing in nuclei 
is still too much model dependent, we need to rely on very accurate quantitative 
theoretical predictions of its properties such as energy gap, superfluidity density 
domain, critical temperature, and other physical quantities associated with the 
various superfiuid states. These quantities may only be obtained from ab initio 
calculations, i.e., microscopic approaches using as input the bare nucleon-nucleon 
interaction, because we are exploring a density domain much wider than the 
saturation region where phcnomcnological interactions such as Skyrme forces 
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are well suited. There is a more basic reason to refrain from using effective 
interactions, that is a double counting of the particle-particle (p-p) correlations 
incorporated in the effective interaction, but also in the gap equation. 

Fortunately, for more than two decades realistic potentials, based on field- 
theoretical approaches, have been supplied to describe the bare nucleon-nucleon 
interaction. The term 'realistic' means that the parameters contained in such po- 
tentials are adjusted to simultaneously reproduce the experimental phase shifts 
of nucleon-nucleon scattering and the binding energies of the lightest nuclei. 

In this chapter the problem of superfluidity in neutron matter is surveyed 
and special emphasis is devoted to new theoretical developments and calcula- 
tions. In the following section the general formalism of pairing in a strongly 
interacting Fermi system is presented. In Sect. || the possible superfluid states of 
nuclcon matter are described within the BCS theory extended to non-zero angu- 
lar momentum. In the last section some aspects of the generalized gap equation 
will be discussed, including the medium polarization effects at very low den- 
sity (Sect. 3.1), the induced interaction approach (Sect. |3^ ), and the role of 
self-energy corrections (Sect. |3.3[ ). 

We mention three previous works for a comprehensive study of superfluidity 
in nuclear matter: the early papers based on the generalized BCS-Bogolyubov 
theory jll],[L2) , which already give a systematic survey of most superfluid states of 
nuclear matter; the second one [HJ, based on the method of the correlated basis 
functions, mainly focussed on the ISo pairing, but containing a wide discussion 
of the important medium correlation effects; and lastly, a recent more general 
overview of pairing in nuclear matter JL4| . 



1.1 Green Function Formalism, Generalized Gap Equation 

In this section we briefly review the main points of the treatment of a superfluid 
Fermi system within the Green function formalism. A detailed account is given 
in various textbooks |l|,|l|,0Jl|,|l|,|§ . 



The principal equations describing a superfluid system are the Gorkov equa- 
tions, that can be considered a generalization of the Dyson equation for a normal 
Fermi system. A diagrammatic representation of the Gorkov equations is shown 
in Fig. |l|(a). They express the relation between normal and anomalous propaga- 
tors G and F, defined by 

G(l,2) = j(T(^)) ) (i a) 

F(l,2) = i(T(!P 1 !P a )>, (lb) 

and the self-energy S and gap function A In a homogeneous system, these four 
quantitites depend only on a four- vector k = (fco, k) and can be written as 2 x 2 
matrices in spin space, for example, 
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Fig. 1. (a) Diagrammatic representation of the Gorkov equations, (b) Equations for 
the self-energy £ and gap function A 



Using the free fermion propagator 
Go(fc) 

and defining 



e(k) 



k a - k 2 /2m + ^ + iOk 
k 2 



2m 



S(k ,k) - fi , 



one can write the system of equations explicitly as 



[k -e(+k)]l A(k) \{G(k) 
At(fe) [ko+e(-k)]l Mrt(fc) 



(3) 



(4) 



(5) 



where 1 denotes the two-dimensional unit matrix. In order to take into account 
at the same time pairing correlations As with spin 5 = and S — 1, one can 
make the ansatz 

+A + iA 1 
-Ao+iAi 



(6) 



and equivalently for F^, whereas the self-energy X and G are diagonal in the 

spin indices. If the ground state is assumed to be time-reversal invariant, the 

gap function has in general the structure of a unitary triplet state ]2l| , p2| |, i.e., 
it fulfills 

A t (fc)A(fc) = Z\(fc) 2 l, (7) 

where by A(k) 2 we denote the determinant of A in spin space. 
The system Eq. (0) can then be inverted with the solution 



G(k) = 
Fl(k) = 



k + e(—k) 

W) 

D(k) ' 



(8a) 
(8b) 
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where 

D{k) = [fc - e(+k)] [fc + e(-k)] - A (k) 2 - A,(k) 2 , (9) 

that expresses the propagators G and F^ in terms of S and A, respectively. 

In order to determine uniquely the four quantities one needs two more equa- 
tions, which relate S and A to the interaction. These equations are displayed in 
Fig. p](b) and read explicitly 



E aa (k) = \J ^^^(fca,fc'/3|T|fca,fc'/3) G^(fc') 



(10a) 



A a/ 3{k) 



d 4 fc' v 
(2tt) 4 ^ 

a',P' 



(ka, -k(3\r\k'a', -k'(3') F a .^ (fc') , (10b) 



where greek letters denote spin indices and T and r are the scattering matrix and 
the irreducible interaction kernel, respectively. Clearly these equations cannot 
be solved in full generality, but one has to recur to some approximation at this 
stage. The simplest, very common, BCS approximation, is to replace T and r 
by the leading term, namely the bare interaction V . In this case the interaction 
is energy independent and the fco integration in Eqs. (10a 10b) can be carried 
out trivially, leading to 



E(k) = 



-I) 2 r 

fe' 



(k,k'\V + 3Vi\k,k') 
-(k,k'\3Vi - V Q \k',k) 



A s (k) = Y l (usv)k>(+k',-k'\Vs\+k,-k) a 
fc 



2 1 (, e 
V = — 1 

2 V E 



2E 



where 



E 2 



2m 



Z{k) 



/J 



(11a) 
(lib) 

(12) 



Together with an equation fixing the chemical potential fi for given density p, 

P = 2J24, (13) 



this is the coupled set of equations that needs to be solved in order to find 
the Hartree-Fock self-energy EnF(k) and the BCS gap function zAbcsW m a 
superfluid system. In contrast to a normal Fermi system, the smooth occupation 
numbers w 2 , instead of the Fermi function 6(kp — |fc|) appear in the HF equation. 



2 BCS Approximation 

In this section we present the solutions of the BCS gap equation 
AMk) = -J2(k\V TS \k')^±, 



k F _ 9 1 



<k) - 
E(k)_ 



(14a) 
(14b) 
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where 



E(k) 2 = e{kf + J2 A Ts{kf , e(k) = e{k) - fi 



(15) 



T,S=0,1 



with fi being the chemical potential and e(k) the single-particle spectrum. Dif- 
ferent realistic nucleon-nucleon potentials V |^,^4 25 2(||27],^8| will be used as 
input. The equations are valid for pure neutron matter (T = 1) and also for 
symmetric nuclear matter (T = 0, 1), having extended the derivation in the pre- 
vious section by including the isospin quantum number T in analogy to the spin 
S. 



2.1 Pairing in Different Partial Waves 



In order to reduce the three-dimensional integral equation ( 14a ) to a set of one- 
dimensional ones, it is advantageous to perform partial wave expansions of the 
potential and the gap function. In this way one arrives at separate equations 
in the different (TSLL 1 ) channels of the interaction, provided an angle-average 
approximation is made by replacing A{k) 2 — > J dk / Att A(k) 2 in Eq. (|l5|). The 
following equations for the partial wave components of the gap function are then 
obtained @Hg|l]|2l||: 



ArsL(k) = — 



y/e(k'f + A(k') 



A2 



-.A 



TSL 



(*') 



(16) 



where 



A{kf = A TSL{kf , 



(17) 



T,S,L 



and with the matrix elements of the bare potential in momentum space 



V[3(k,k')= I drr i j L ,(Vr)V 1 



l3(r)j L (kr) 



(18) 



It should be noted that the different equations are still coupled due to the fact 
that the total gap appearing in the denominator on the r.h.s. of Eq. (16) is the 
r.m.s. value of the gaps in the different partial waves. The gap equation allows in 
principle the coexistence of pairing correlations with different quantum numbers 
(TS), even though the different (TS) channels are not mixed by the interaction. 
In practice, however, so far no such mixed solutions of the gap equation have 
been found: even if at a given density two or more uncoupled solutions exist, the 
strong nonlinear character of the gap equation prohibits a coupled solution. (In 
finite nuclei such mixed solutions seem to exist under certain conditions ]34|). 
This means that in practice this (T5*)-coupling can be neglected and that at a 
given density the solution of the uncoupled gap equation with the largest gap is 
the energetically favoured one. 

The only case when it is clearly necessary to keep the coupled equations is 
the mixing of partial waves due to the tensor potential. In this case the gap 
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equation can be written in matrix form (for given S = 1 , T, L) 



A L 
Al, 



(k) 



1 f°°dk'k' 2 1 ( Vll Vll + 2 

E(k') \V L+2 . L V L+ 2,L+2 



(k,k') 



A L 

A L+2 



(k') 

(19) 



with 



E(k) 2 = [e(fc) - M ] 2 + A L {k) 2 + A L+2 {k) 2 . (20) 

This is relevant equation for the 3 SD\ (T = 0) and 3 PF 2 (T = 1) channels, for 
example. 

Let us finally mention that usually, apart from the 3 <SDi channel, the two 
equations (14a) and (14b) can be decoupled by setting /x = e(kp). The reason 
is the small value of the ratio Aj /i, so that a Fermi surface is still quite well 
defined. 

We come now to the presentation of the results that are obtained by solv- 
ing the previous equations numerically, using a kinetic energy spectrum e(fc) = 
k 2 /2m for the moment. In practice, one finds in pure neutron matter (T = 1) gaps 
only in the % P,pH35|, |3^,p7||38| a nd 3 PF 2 [[|j2|,|^,|l|j3|j3| partial waves. 
They are reported in Figs. |2| and |3l respectively. It can be observed that the 
maximum pairing gap is about 3 MeV in the channel and of the order of 
1 MeV in the 3 PF 2 wave. It is remarkable that solutions obtained with different 
nuclcon-nucleon potentials are nearly indistinguishable in the case, whereas 
in the 3 PF 2 wave such good agreement can only be observed up to kp m 2 fm _1 
(with the exception of the Argonne Vxa potential that is not very well fitted 
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Fig. 2. So gap evaluated in BCS approximation with free single-particle spectrum and 
different potentials 



Superfluidity in Neutron Star Matter 7 

P [fm 3 ] 



0.0 0.1 0.2 0.4 0.6 0.8 1.0 1.5 




k F [fm" 1 ] 

Fig. 3. 3 Pi<2 gap evaluated in BCS approximation with free single-particle spectrum 
and different potentials 



to the phase shifts), from where on the predictions start to diverge from each 
other. The reason 13T 



J8j is the fact that the various potentials are constrained 
by the phase shifts only up to scattering energies E\ a h of about 350 MeV, which 
roughly corresponds to a Fermi momentum of kF ~ \J "i-E-iab/2 ~ 2fm~ 1 . Thus 
even on the BCS level the gap in the z PFi channel at a neutron density higher 
than w 0.3 fm~ 3 is at the moment not known. Apart from that it is clear that 
the BCS approximation is not reliable at the very large densities for which a 
gap is predicted in Fig. |3|. However, nobody has so far attempted to include 
polarization effects in this channel. 

Let us mention here for completeness that in symmetric nuclear matter one 
finds very strong pairing of the order of 10 MeV in the Z SD\{T = 0) channel, 
reminiscent of the deuteron bound state @|T]J2||^,|^I^2[||,||,||J4§1 , and 
also a gap of the order of lMeV in thc^D 2 wave This, however, is 

probably not very relevant for neutron star physics, since a prerequisite for this 
T = neutron-proton pairing to take place is the existence of (nearly) isospin 
symmetric nuclear matter, the pairing correlations being rapidly destroyed by 
increasing asymmetry [ ffSp^Jsof . We will therefore not discuss this type of pairing 
further on. 



2.2 Pairing Gaps in Neutron Star Matter 



Let us now come to the T = 1 gaps that can be expected in isospin asymmetric 
(beta-stable and charge neutral) neutron star matter. On the BCS level, the only 
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Fig. 4. The different T = 1 gaps in neutron star matter as a function of total nucleonic 
density 



influence of isospin asymmetry on these gaps is via the neutron single-particle 
energy e(k) appearing in the gap equation. Calulations have been performed 
using e(k) determined in the BHF approximation extended to asymmetric nu- 
clear matter ||M||,||,|3[||g|§, and typical results are shown in Fig. |. 
One observes that Sq pairing can take place independently in the neutron and 
in the proton component of the matter. If plotted as a function of total baryon 
density, the neutron pairing occurs naturally at lower density and with a larger 
amplitude than the proton pairing, because at the higher density the proton 
effective mass is smaller and the pairing therefore more reduced. The same is 
true for the 3 PF 2 pairing in the neutron component, which is strongly reduced 
with respect to the calculation with a free spectrum shown in Fig. || above. We 
stress that the results displayed in Fig. ^ can only be qualitative, because they 
clearly depend on the details (in particular the proton fraction) of the equation 
of state that is used. The results shown were obtained with a BHF EOS based 
on the Argonne V14 potential and involving n, p, e, /1 components p2]. 



3 Beyond BCS 

In the previous section we have presented many results that were all obtained 
within the BCS approximation. However, as has been explained in the intro- 
duction, this approximation amounts to a mean-field approach, equivalent to 
and consistent with the Hartree-Fock approximation in a normal Fermi system. 
More precisely, the BCS approximation neglects completely any contribution 
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beyond the bare potential to the interaction kernel T appearing in the general 
gap equation (10b). 

Going consistently beyond the BCS approximation is however a very difficult 
task and has been only partially achieved so far. We will review in the following 
sections some aspects of these extensions. We begin with a discussion of the 
situation at extremely low density, where certain analytical results are known. 
Following that, the general framework at more relevant densities will be set up, 
and we will briefly present the results that have been obtained so far by various 
authors. Finally, in the last section, we will focus on a certain part of the problem 
that has recently been tackled, namely the treatment of the energy dependence 
of the self-energy that appears in the gap equation. 



3.1 Low Density 

In order to derive an exact analytical result for the pairing gap including polar- 
ization effects that is valid at very low density (more precisely, for kp « l/\a\, 
where a is the relevant scattering length), we begin again with the BCS gap 
equation, 

A k = -Y,V k k'^ r A k , , E k = ^f(e k ~e F y + Al. (21) 
fc' k ' 

It is then useful J3(|^4|,^7| to introduce a modified interaction T that is given by 
the solution of the integral equation 

Tkk' = Vkk' — ^^T kk nF k "Vk"k' , (22) 
fc" 

where F k is for the moment an arbitrary function. (We have used the symbol 
T, although in general this quantity is not to be identified with the scattering 
matrix). Making use of this equation, the gap equation is transformed into 

A k = -J2 T kk' -(^r; ~ Ffe ') Ak ' ■ ( 23 ) 

Of particular interest is now the choice F k = sgn(fc — kp)/2E k , which leads to 
the set of equations 

A k = -2 ]T T kk ,^-A k ,, (24) 
fc'<fc F k ' 

T kk ' — V kk > — y^Tfcfc" — — — —V k » k > . (25) 

Therefore, in the limit A/ep — > that is approached with vanishing density, 
T does become identical to the free scattering matrix, because E k — > \e k — ep\ 
in this situation. At the same time, in the gap equation ( p4| ) the interaction T 
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is now cut off at k' = kp (at the cost of introducing a factor 2), so that with 
vanishing density it is ultimately sufficient to use the low-energy result |5S[ ] for 
the T-matrix: 

T kk , - r 00 = ^ , (26) 
m 

where a nn = — 18.8fm is the neutron-neutron scattering length. This yields 
finally the gap equation 



1 = - 



Ak F a nn 



dx ■ 



y/(l-x*)* + (A/e F y 



which in the limit A/e F -> is solved by |7jj59|j6C|,|6lJ 

A(k F ) k ^° z\ (fc F ) = ^M:exp 
e z 2m 



2k F a r , 



(27) 



(28) 



This is the universal asymptotic result for the BCS pairing gap in a low-density 
Fermi system with negative scattering length. Unfortunately its validity is lim- 
ited to the region k F « l/\a nn \ « 0.05 fm _1 , far below the densities of interest 
for neutron star physics or even pairing in finite nuclei. 

Going now beyond the BCS approximation, in the low-density limit one 
should take into account the corrections to the interaction kernel that are of 
leading order in density. Diagramatically these are the polarization diagrams 
of first order (i.e., comprising one polarization "bubble") that are displayed in 
Fig. H The interaction appearing in these diagrams is in the present case the 
free scattering matrix T. It can be shown |5{| that in the low-density limit it 
is again sufficient to neglect the momentum dependence of the T-matrix, as in 
Eq. (E6J) . One obtains then for the lowest-order polarization interaction 



Wi Sn (k,k') = -^M — 
5oV ' ; 8tt 2kk' 



k+k' 



dqqU{q) 



\k-k'\ 



(29) 



where 



mk F 

1%) = — 

TV* 



1 1 - X 2 



4.c 



In 



1 + x 



1 - x 



* = w r < 30 ' 



! v 
4 — 





(a) 



(b) 



Fig. 5. (a) Bare potential and (b) first-order (direct and exchange) polarization dia- 
grams contributing to the interaction kernel in the low-density limit 
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is the static Lindhard function p 
the BCS interaction kernel Eq. ( 



Thus the lowest-order polarization modifies 
J) by a (repulsive) term proportional to kp, 



whereas any other polarization diagram contributes only in higher order of Uf- 
We can use this result and insert it into the previously obtained approxima- 
tion in the BCS case: 



kp 
2m 



exp 



tt/2 



where 



2tt 
mkp 



2kF Aqq 2 

— n(g) = — (1 



2kj 



3tt 



kpa n 



2 In 2) « 0.506 



(31) 



(32) 



accounts for the polarization effects to first order. Expanding now the argument 
of the exponential up to second order in k, one obtains for the ratio relative to 
the BCS value, Eq. (§8|), 



A{k F ) 
A (k F ) 



exp 



1 



-~c[1-ck + G{k 2 )] 

(1-ck) 



(4e)V3 

fe F ^0 1 

"~~ * (4e)V3 ■ 



(33) 

(34) 
(35) 



Let us stress that the above "derivation" of this result can only be considered 
heuristic. A rigorous proof was given originally in Ref. [ p9[ . 

Therefore, one arrives at the striking conclusion that in the low-density limit 
the polarization corrections suppress the BCS gap by a factor (4e) -1 / 3 w 0.45, 
independent of the strength of the interaction a nn . This is quite surprising, since 
the polarization interaction clearly vanishes with vanishing density; its effect on 
the pairing gap, however, does not. The reason is the nonanalytical dependence 
of the gap on the interaction strength, as expressed by Eq. (|28|) . 

All this means that the BCS approximation cannot even be trusted at very 
low density, and that one can in general expect quite strong modifications due 
to polarization effects at higher density as well. 

This low-density behaviour of neutron pairing could be met in the study of 
exotic nuclei with a long density tail or in nuclei embedded in a neutron matter 
environment, as occurring in the neutron star crust. 



3.2 General Polarization Effects 

To go beyond the simple low-density approximations derived before requires 
considerable effort and has in fact so far not been accomplished in a satisfactory 
manner, so that ultimate results cannot be presented here. The reason is that 
many effects that could be neglected in the low-density limit become important 
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now, and that on the other hand the pairing gap is extremely sensitive to even 
slight changes of the interaction kernel. 

First, outside the low-density region Uf « l/|a nn |, the interaction kernel 
has to be extended beyond the lowest-order polarization diagrams. This means 
summing up polarization diagrams of all orders in the particle-particle channel, 
but also including them in the particle-hole channel, replacing the T-matrix by 
the general particle-hole interaction F. In this way a self-consistent scheme is 
established that is depicted in Fig. [| It requires as input the Brueckner G-matrix 




Fig. 6. Determination of the interaction kernel F in the gap equation (a): Polarization 
diagrams appear in the particle-particle channel (b) as well as in the particle-hole 
channel (c). The leading diagrams in these channels are the bare potential V (dashed 
line) and the G-matrix (double-dashed line), respectively 

and yields ideally the interactions in the particle-particle as well as particle-hole 
channel, J 1 and F, respectively. It is clear that in practice an exact solution is 
impossible, but that usually strong approximations have to be performed that 
cast a doubt on the reliability of the results that are obtained. We will later 
discuss this point in some more detail. 

Second, related to the previous item, the energy dependence of the full gap 
equation [see Eq. ( |l0b|) 1 needs to be taken into account. This is obvious, since the 
interaction kernel becomes now a complex, energy-dependent quantity. To our 
knowledge, this problem has so far not been studied in detail in the literature. It 
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is therefore not known how far the gap could be changed by this more elaborate 
treatment of the equations. 

Third, and in connection with the two previous points, the choice of a par- 
ticular interaction kernel requires also the choice of a compatible self-energy 
appearing in the gap equation. This will be explained in more detail in the 
following section, as it has recently been addressed in the literature. 

It should be clear by now that the influence of medium effects on pairing 
constitutes an extremely difficult problem. Consequently the results that can be 
found in the literature p^| , |62| , |63|j6^ , |65| , |66[ addressing this task in certain approx- 
imations agree only on the fact that generally a strong reduction with respect 
to the BCS gap is obtained. A collection of these results is displayed in Fig. ^. 
It can be seen that the precise amount and density dependence of the suppres- 
sion vary substantially between the different approaches and must be considered 
unknown for the time being. 

The most advanced description of medium polarization effects is based on 
the Babu-Brown induced interaction model J6^ , |68|j6^ , |70| , [7l| , [72[ . The microscopic 
derivation of the effective interaction starts from the following physical idea: The 
particle-hole (p-h) interaction can be considered as made of a direct component 
containing the short-range correlations and an induced component due to the 
exchange of the collective excitations of the medium. 

Let us consider a homogeneous system of fermions interacting via an instan- 
taneous potential V , which is also translationally and rotationally invariant. Col- 
lective excitations are described by the ring series, which can easily be summed 
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Fig. 7. The So gap in pure neutron matter predicted in several publications taking 
account of polarization effects. The curve in the background shows the BCS result 
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up p8[ , and the p-h interaction can be written as (disregarding for the moment 
spin degrees of freedom) 

In this simple case the interaction itself plays the role of the direct term and the 
sum of the ring series that of the induced term. 

In the nuclear case the presence of the hard core imposes the bare interaction 
V to be renormalized in order to incorporate the short-range correlations. This 
goal is reached by introducing the G-matrix, which sums particle-particle (p-p) 
ladder diagrams to all orders, and the diagrammatic expansion can be recast 
just replacing V by G. But now the new ring series cannot be summed up any 
longer, mainly since the G-matrix is nonlocal. An averaging procedure has been 
devised to bring the G-matrix into a local form G{q) |6(| . Then, in analogy to 
Eq. (|36|), the p-h interaction F is given by 

In a simplified version of the theory the direct term now coincides with the G- 
matrix and the induced term is the approximate sum of the renormalized ring 
diagrams. In Fig. || the direct term is represented by the first diagram on the 
r.h.s. of the series (c). The next diagrams form the ring (or bubble) series. 

But, since the RPA series with the G-matrix produces a too strong polariz- 
ability of nuclear matter, in the Babu-Brown approach it has been proposed to 
include in the RPA series the full p-h interaction itself, since the particle (hole) 
coupling vertex with the p-h bubble can indeed be identified with the irreducible 
p-h interaction. The series (c) of Fig. || is, in fact, the final result of such a proce- 
dure, where the wiggles represent the effective p-h interaction F on either side. 
If we denote by Fd the direct interaction (G-matrix in our case) and by Fi the 
induced interaction, the effective p-h interaction can be written in the form (at 
the Landau limit, in which p\ and p 2 are restricted to the Fermi surface) 

F{px,pi) =F d ( Pl ,p 2 ) + F i {p 1 ,p 2 ;F), (38) 

where Fi is, as said before, the RPA series of Fig. |^(c) with the G-matrix replaced 
by F itself. The previous equation clearly entails a self-consistent procedure to 
determine the interaction F{p\,p2). Referring to |6(| for details, one may reduce 
Eq. ( |38| ) to a numerically tractable form 

which corresponds to replacing the G-matrix by F in the induced term. 

Once the irreducible p-h interaction F has been determined, one can con- 
struct the irreducible p-p interaction r by performing the tranformation of the 
matrix elements of the interaction from the p-h to the p-p channel. However, 
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this is not enough, since in the p-p channel a set of additional diagrams must 
be added arising from p-h diagrams which are reducible in that representation. 
After including these terms, the interaction contains both direct and exchange 
terms, which guarantees antisymmetry and Landau sum rules and, in addition, 
it should simultaneously make the nuclear matter Landau parameter Fq less 
negative so that the stability condition is satisfied. The first contributions to 
the p-p interaction are depicted in line (b) of Fig. [| We stress once more that 
they describe the influence of the medium polarization on the nucleon-nucleon 
interaction within the induced interaction model of Babu-Brown. 

We discuss now the effects of medium polarization on the superfluidity of 
neutron stars in the channel ^So- First of all, the irreducible p-p interaction 
to be used in the pairing problem must not include any ladder sum already 
included in the gap equation, and therefore the first term in line (b) of Fig. || is 
the bare neutron-neutron interaction V. The next terms include the irreducible 
p-h interaction in the vertices of the p-h bubble. If their momentum dependence 
is neglected, these vertices can be identified with the Landau parameters (6j| . 

In neutron matter the polarization of the medium is due to density fluctua- 
tions and spin-density fluctuations given by 



respectively. Solving the Babu-Brown self-consistent equation, Eq. (|39|), with 
the G-matrix as direct interaction and the renormalized RPA series as induced 
interaction, one determines the p-h interaction and eventually, after the Landau 
angle expansion, the lowest-order Landau parameters Fq, related to the nuclear 
compression modulus, and Go, related to the spin waves. (In the following the 
interaction will be expressed in terms of these two Landau parameters for sim- 
plicity). Then the effective interaction is calculated including the diagrams of 
Fig. |6|(b). The pairing interaction in the 1 Sq channel is then given by 



This equation shows that the medium screening effect is determined by the 
competition between the attractive term induced by density fluctuations and 
the repulsive term induced by spin-density fluctuations (first and second term 
in the bracket, respectively). 

As depicted in Fig. ||, the effect of the medium polarization is an overall 
suppression of the gap due to the prevalence of the spin-density fluctuations 
over the density fluctuations. A similar effect is found in a less crude calculation 
|p6| , as shown in Fig. |7[ where the peak value is shifted to higher density. Such a 
suppression is common to all calculations existing in the literature even if, as to 
its magnitude, the different predictions do not agree with each other, as shown 



5pk = &Pk] + Spki , 
Spk = Spk^ - Spki , 



(40a) 
(40b) 




(41) 



in Fig. 0. 
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Fig. 8. So pairing gap in neutron matter as a function of the Fermi momentum kp. 
The full curve corresponds to using the bare V14 potential as direct interaction. The 
symbols show the effect of the medium polarization described in terms of the Landau 
parameters, according to Eq. mW) 



3.3 Self-energy Effects 

Dynamical effects of the interaction on the gap function have been completely left 
aside from the discussion on the medium polarization. So far no solution of the 
gap equation has been attempted considering the irreducible interaction block 
as an energy-dependent quantity. On the other hand the energy dependence in 
the self-energy can affect deeply the magnitude of the energy gap in a strongly 
correlated Fermi system such as nucleon matter ]73| , |74| |. 

To discuss self-energy effects we come back to the generalized gap equation 
presented in section [hi]. Let us rewrite Eq. (10b) in the following form (for the 

channel) 

A , . f d 3 k' fdu' . n A k ,(u') 

4(w) ^i^i^ %(w ' w) ^M' (42) 



with [cf. Eqs. (pah and 



- D k {io) = [G k {-Lo)Gt{u)]- 1 = Gi\«>)G?{-w 



+ (43) 



The functions G k {uo) and Gt(u>) are the nucleon propagators of neutron matter 
in the normal state and in the superfluid state, respectively. The w-symmetry in 
the two propagators is to be traced to the time-reversal invariance of the Cooper 
pairs. The effective interaction r is the block of all irreducible diagrams of the 
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interaction. Since we want to focus only on the self-energy effects, we assume the 
interaction to be the bare interaction Vk,k' , as in BCS. Then the pairing gap does 
not depend on the energy (static limit), i.e., A k {oS) = Ak, and the w-integration 
can be performed in Eq. (|4^), once the self-energy has been determined. A general 
discussion of the analytic w-integration of the gap equation has been is given in 
Ref. [jfij. Here we follow a simplified treatment based on the fact that, at each 
momentum fc, the main contribution to the w-integration comes from the pole 
of Gk{uj). This latter is the solution of the implicit equation 

k 2 /2m + Skiuifs) — \i ■ 



(44) 

Expanding the self-energy around the pole ui k amounts to expanding G _1 itself, 
which yields 

'-ilj^-i' <45) 

where the prefactor on the r.h.s. is the inverse of the quasiparticle strength Z k 
that will be discussed later. Then the energy dependence of the kernel (w) 
takes the simple form 



zl 



- W 2 -(Z fe zA fc ) 2 ' 

and the w-integration can be performed in the usual way, leading to 

d 3 fc' ,„ „ „ , Ay 



A k = 



(2tt): 



2^fu 



(46) 



(47) 



Comparing with the BCS result, Eq. (14a), the new gap equation contains the 
quasiparticle strength Z to the second power. Since Z is significantly less than 
one in a strongly correlated Fermi system, a substantial suppression of the energy 
gap is to be expected. 

Moreover Z deviates from unity only in a narrow region around the Fermi 
surface and hence it is not a severe approximation to restrict the w-integration 
to only the pole part at the Fermi energy. Expanding then £k(u>) around the 
Fermi surface (u> = and k = hp), Eq. ( pfjj ) is easily solved and we get ui k ~ 
(fc 2 — kp)/2m*, where m* is the effective mass at kp [see Eq. (f53|)]. The gap 
equation becomes 



A, 



7 2 



d 3 fc' 



Vkk'Aii 



(48) 



(2tt) 3 2 N /[(fc' 2 - k' F )/2m*} 2 + A{, ' 

where Zp is the quasiparticle strength at the Fermi surface. 

As is well known the pairing modifies the chemical potential which is calcu- 
lated self-consistently with the gap equation from the closure equation for the 
average density of neutrons 



P 



Z F 



d A k 

d 3 k 



27ri 



(2tt) 3 



(fc 2 



kp)/2m* 



^[(fc 2 - k; F )/2m*f + Z\ 



A 2 



(49) 



(50) 
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The latter approximation is sufficient to investigate the self-energy effects. 

Before we present the predictions based on Eq. ( [48| ) , we discuss some proper- 
ties of the self-energy £k(&) of neutron matter. In the Brueckncr approach |75| 
the perturbative expansion of S can be recast according to the number of hole 
lines as follows 

^H = 4 1) h + 4 2) m + ... . (5i) 



G 




Fig. 9. Hole-line expansion of the self-energy 

The on-shell values of represent the Brueckner-Hartree-Fock (BHF) mean 
field (first diagram in Fig. |); the ones of S^ 2 ' represent the so-called rearrange- 
ment term (second diagram in Fig. which is the largest contribution due 
to ground-state correlations. The off-shell values of the self-energy are required 
to solve the generalized gap equation. Fig. 10 displays a typical result for the 
off-shell neutron self-energy Ek(uj) calculated up to the second order of the hole- 
line expansion. The calculations are based on Brueckner theory adopting the 
continuous choice as auxiliary potential |76| . 

Since we are interested in the behaviour of the self-energy around the Fermi 
surface (k — kp and u> ~ 0), we may use the expansion 



S k (u) «^ fcF (0) + 



8£_ 

duo 



dS 
~dk 



(k - k F ) 



(52) 



0.8 Pm 



■l 



k F = 0.8 Frrf 1 
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co - co F [MeV) 



-50 50 
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i ■ ■ 1 1 i 
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-50 50 
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Fig. 10. Off-shell self-energy in neutron matter at kF = 0.8 fm - . Solid curves: BHF 
approximation Dashed curves: rearrangement contribution X 1 ' 2 '. The Argonne 

V14 potential |24| has been used in the Brueckner calculations 
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In this case the quasiparticle energy takes the simple form 
k 2 ~k 2 F l + (m/k F )(d£/dk)\ F k 2 - 



k 2 F 



2m 



1 - (dE/dco)\j 



2m* 



(53) 



In the previous equation we have introduced the effective mass m*/m as the 
product of the e-mass m e and the fc-mass m kl which are defined respectively as 
follows |73] 



m e 
m 

m k 
m 



1 



duo 



1 



_^ m dS 
kp dk 



(54a) 
(54b) 



The partial derivatives are evaluated at the Fermi surface. The fc-mass is related 
to the non-locality of the mean field and, in the static limit (w = u>f), coincides 
with the effective mass. This quantity is of great interest in heavy-ion collision 
physics, since the transverse flows are very sensitive to the momentum depen- 
dence of the mean field. The e-mass is related to the quasi-particle strength. 
This latter gives the discontinuity of the neutron momentum distribution at the 
Fermi surface, and measures the amount of correlations included in the adopted 
approximation. 

From Ek(w) the effective masses are extracted according to Eqs. ( |54a| ) and 
( ^4b| ). They are depicted in Fig. |ll|, where the full calculation is compared to 



k F = 0.65 Pm" 1 
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Fig. 11. Effective masses in neutron matter for three densities: effective mass m* (up- 
per panels) and e-mass (lower panels). Dashed lines correspond to E^ and solid lines 
to E^ + E^ 
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that including only the BHF self-energy. We may distinguish two momentum 
intervals: at k m kp the momentum dependence of the effective mass m* is 
characterized by a bump, whose peak value exceeds the value of the bare mass; 
far above kp the bare mass limit is approached. The contribution from the 
rearrangement term exhibits a pronounced enhancement in the vicinity of the 
Fermi energy, which is to be traced back to the high probability amplitude for 
p-h excitations near ep frf|. At high momenta this contribution vanishes. One 
should take into account that in this range of kp the neutron density is quite 
small (p = 0.074fm~ 3 at the maximum kp = 1.3 fm -1 ). This behaviour of the 
effective mass m* is mostly due to the e-mass, as shown in the lower panels of 
Fig. [ll]. In all panels of Fig. O is also reported for comparison the effective mass 
in the BHF limit (only Z^ 1 ) included), which exibits a much less pronounced 
bump at the Fermi energy. 

With the off-shell values of the self-energy discussed above as input, the 
gap equation has been solved in the form of Eq. (48), coupled with Eq. ( |50| ) 
[ frr| . This is a quite satisfactory approximation, especially in view of studying 
the self-energy effects on pairing. The Argonne V14 potential has been adopted 
as pairing interaction, which is consistent with the self-energy data for which 
the same force has been used. The results are reported in Fig. [l^ for a set of 
different kp -values. The diamonds connected by a solid line represent the solution 
of Eq. p8| ) replacing the effective mass by the free one, m* — m, together with 
Z = 1 . This is very close to the prediction obtained from the BHF approximation 




0.00 0.25 0.50 0.75 1.00 1.25 1.50 
k F [frn -1 ] 

Fig. 12. Energy gap in different approximations for the self-energy: free s.p. spectrum 
(solid line); effective mass with Zf = 1 (upper dashed line); Z from X 1 ' 1 ' (long dashed 
line); Z from E m + £< 2 > (lo wer dashed line) 



Superfluidity in Neutron Star Matter 21 



for the effective mass m* , but still keeping Z = 1 (fancy crosses connected by 
the dashed line). This similarity stems from the fact that at the Fermi surface 
m*/m from BHF is close to one, as shown in Fig. [n]. The self-energy effects 
are estimated in two approximations. In the first one m* and the Z-factor are 
calculated from a BHF code. In the considered density domain the Z-factor 
ranges around 0.9. Despite its moderate reduction a dramatic suppression of the 
gap is obtained, as shown by the long dashed line in Fig. [l^. It is due to the 
exponential dependence of the gap on all quantities. Still a further but more 
moderate reduction is obtained when the rearrangement term is included in the 
second approximation. The smaller Z-factor (Z « 0.83 at &f = 0.8fm _1 ) is to a 
certain extent counterbalanced by an increase of the effective mass (to* /to w 1.2 
at the same ftp). 

From the previous discussion we may conclude that including self-energy 
effects is an important step forward in understanding the pairing in nucleon 
matter, but one should also include, on an equal footing, vertex corrections. 

4 Conclusions 

This chapter addressed the present status of the theoretical progress on the 
superfluidity of neutron matter and the microscopic calculations of nucleonic 
pairing gaps relevant for the conditions that are encountered in the interior of 
a neutron star. We have seen that most of the results that can be found in the 
literature, are obtained within the BCS approximation. 

Unfortunately, as has also been pointed out, this approximation cannot be 
considered reliable in any region of density, since neutron matter is a strongly 
correlated Fermi system: The same nucleons that participate to the pairing cou- 
pling also participate to screen the pairs. This requires necessarily to improve 
the interaction kernel by in-medium polarization corrections that, at the present 
time, can be treated with confidence only at very low density, where the ring 
series can be truncated at the order of the one-bubble term. 

The induced interaction approach is a promising candidate to describe these 
effects, but its predictions are still affected by too severe approximations. The 
existing calculations addressing this problem point to a substantial reduction of 
pairing in the (T = 1) 1 Sq channel with respect to the BCS results. In the 3 Pi*2 
channel, no estimate of polarization effects has been made so far, at all. 

Also self-energy effects, when properly included in the gap equation, strongly 
influence the pairing mechanism. We saw in fact that the gap is reduced with 
the square of the quasi-particle strength Z. 

In any case it remains the problem of a complete solution of the generalized 
gap equation, taking into account simultaneously screening and self-energy cor- 
rections, which moreover have to be treated on the same footing. This latter goal 
remains a considerable theoretical challenge for the future. 

Due to lack of space and time we have not discussed more speculative subjects 
like pairing in isospin asymmetric matter, relativistic effects on pairing, hyperon 
pairing, etc., which might also have relevance for neutron star physics. 
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